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ABSTRACT 

We find evidence that the mass Mbh of central supermassive black holes (SMBHs) 
correlates with the velocity dispersion ctqc of globular cluster systems of their host 
galaxies. This extends the well-known Mbh — <?sph relation between black hole mass 
and velocity dispersion of the host spheroidal component. We compile published mea- 
surements of both Mbh and ctqc for a sample of 13 systems and find the relation 
log(M BH ) = a + /31og(cr GC /200) with a = 8.63 ± 0.09 and j3 = 3.76 ± 0.52. We also 
consider blue (metal-poor) and red (metal-rich) globular clusters sub-populations sep- 
arately and obtain a surprisingly tight correlation using only the velocity dispersion 
Oqc of the red clusters with a = 8.73 ± 0.09 and j3 = 3.84 ± 0.52 and an intrinsic 
scatter Eq = 0.22 dex compared to Eq = 0.27 dex for the Mbh — <?sph of our sample. We 
use this Mbh — o-qc relation to estimate the central black hole mass in five galaxies 
for which (Jqq is measured. 

Key words: black hole physics - galaxies: evolution - galaxies: nuclei - galaxies: 
star clusters: general - galaxies: fundamental parameters - globular clusters: general 



1 INTRODUCTION 



It has been known for more than a decade that supermassive 
black holes (SMBHs) reside at the centre of many galaxies 
ijKormendv fc R ichstonc 199$|) and that their mass corre- 
lates with various properties of the ir host galaxies, for ex - 
ample with the spheroid luminosity dMagorrian et al1ll99ct). 



12009] , IVolonteri fc Nataraiarj|2009l ). Furthermore, the phys- 
ical scale at which the feedback is effective is still 
not clear. Studies emphasising that black hole mass 



than the spheroid (Fcrrarcsc [20021: 


Booth & Schavd l20ld: 


IVolonteri, Nataraian, Gultekmll2011 


) seem to indicate that 



the s phero id velocity dispersion a sv h (|Ferr a rose & Merritt 
20001 and iGebhardt et alj l200d . see als o iGultekin et al 



2009 , hereafter G09) , the sp heroid mass (|Magorrian et al 
19981 : iMarconi fc Hunt 2003) an d the kinetic energy of ran- 



dom motions (|Feoli fc Manciml 120091 ). Recently, extension 
of these c o rrelations to ba rre d galaxies and pseud o-bulges 



|Hul 120081 : iGrahaml l2008al lbl: iGraham et al 



201 ll) as well 



as ga l axies ho s ting nuclear star clusters ( Ferrarese et al.l 
l200d : IGrahaml I2OI2I ) has also been proposed. In addi- 
tion, a new correlation between the mass of SMBHs and 
the observed total number of globu lar clusters of their 
host galaxies has been discovered jBurkert fc Tremainel 
20ld, hereafter BT10. iHarris fc Harris! l201ll , see also 



Snyder. Hopkins, fc Herncmistl |2011| ). The physics under 



lying these observations is not fully understood, but it 
is widely believed tha t feedback proces ses are responsible 
for such correlations (|Silk fc Reesl 1 19981 ) . It is not known 
whether these relations evolve with time or are primordial 
and have been set at the time of formation of the SMBHs and 
galax i es jDi Matteo. Springel. fc Hernquistll2005l. ICox et al 



l200fj , iHopkins et al 1 120071 . , I Hopkins, Murray, fc Thompson 



feedback processes operate on large scales even though 
iKormendv fc Benderl (|201lh argue that there is no direct 
relation between DM halos and central black holes. Bet- 
ter and more data, extending to higher redshifts, shall shed 
more light on these questions. 

In this paper we provide evidence that the classical 
Mbh — <Js P h relation between black hole mass and spheroid 
velocity dispersion is also working with the velocity disper- 
sion of globular cluster (GC) systems. We use the published 
data for the velocity dispersion <tgc of the GC systems and 
the estimated black hole mass in 13 galaxies. We then study 
separately the two sub-populations : blue, metal-poor GCs 
and red, metal-rich and younger GCs which are generally 
closer to the centre of the galaxy. We show that the correla- 
tion degrades from the red to the blue GCs sub-population. 

This paper is organised as follows. In Section [2] we 
present the data. The correlation between the mass of 
SMBHs and the velocity dispersions of the GC systems of 
their host galaxies is presented in Section [3] The red and 
blue sub-populations are studied separately in Section [3] In 
Section[5] we estimate the black hole mass for 5 galaxies and 
conclude in Section [6] 
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2 THE DATA 

We have found 19 galaxies for which measurements of the 
line-of-sight velocity dispersion of the GC system are known. 
Thirteen of these are accompanied by the estimate of the 
mass of their central black hole. These galaxies are pre- 
sented in Table[T] To our knowledge, the mass of the central 
black holes have not yet been determined for the six other 
galaxies for which we ha ve the velocity dispersion of the 
GC s ystem : NGC 1407 (|Romanowskv et al.ll2009l) NGC 
3923 (iNorris et al.l 120121), NGC 4 494 jFoster et all [ioTlh . 
NGC 4636 llLee et alj|20ich. LMC [Freeman 1 19931 ) and M33 
l|Schommer et al.lll99ll ; IChandar et al.ll2002T ). 



2.1 Velocity dispersion of globular cluster systems 

Galaxies for which we found measurements of ctgc as well 
as o'gc an< l °Gc e are listed in Table [T] along with the cor- 
responding references. Whenever available, we use the value 
of ctgc corrected for the global rotation of the GC system. 
Comments on individual systems are given below : 



• NGC 1399 - We use the results of iRichtler et al. I l2004l 
which give c tgc for all GCs and for both red and blue sub- 
populations. [Schuberth~eTjiL] (|2010t ) give detailed results for 
different sub-samples and dif ferent distance s whic h are in 
general not too different from IRichtler et"aL I (|2004h . 

NGC 3031 (M81) - Re s ults h a ve been obtain e d sep - 



arat ely by iPerelmuter et all (119951 ) . ISchroder et all (2002) 
an d iNantais fc Huchra (I201C ) . We use the values given 
bv lNantais fc Huchral ()201Ch which are consistent with the 
other groups. Measure ment uncertain ties in ctgc, o"gc and 
CT G ^ e were provided bv lNantais! (|2012h . 

• NGC 4594 (M 104) - The value of ctqc is taken from 
iBridges et all (|2007t ). They also measured ctJJc and ct G c c but 
do not give the corresponding error bars. Consequently, we 
do not include M104 in our best-fit estimate when treating 
red and bl ue GC s sepa rately. 

• MW - IZinnl (1996) gives ctgc for three groups of GCs 
lying approximately around the direction of the galactic cen- 
tre. The corresponding line-of-sights are therefore nearly 
parallels and we can treat the measured radial velocities 
as projection on a single line-of-sight similarly to what is 
done for external galaxies. The three groups are metal-poor 
(blue), red clusters lying between 2.7 to 6 kpc from the galac- 
tic centre, and very metal-rich disc clusters. We choose the 
dispersion of the red clusters (uqq = 61 ± 10 km s _1 ) to 
compare with the dispersion of red GCs in others galax- 
ies. W e note that a similar value was found by lHarrisI 
(1999) for metal-rich clusters lying in the range 4 — 9 kpc. 
Concerning the metal-poor cluste rs, we retain the value 
ct g ^ c = 120 ± 14 km s _1 given by lHarrisI (|l999T ). However, 
we point out that comparison with external galaxies is not 
straightforward in this case since we can no longer consider 
the projection of cluster velocities on a single line-of-sight. 



2.2 Mass of central black holes 

The mass of central black holes in 9 out of the 13 galaxies in 
our sample, namely NGC 224, NGC 1399, NGC 3031, NGC 
3379, NGC 4486, NGC 4594, NGC 4649, NGC 5128, NGC 
7457 and MW have been studied and analysed by G09. To 
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Figure 1. Mass Mbh as a function of ctgc f° r the 12 galaxies 
in our sample. Circles indicate galaxies included in the fitting 
procedure. NGC 4636, for which only an upper limit of Mbh is 
available, is shown as a square. The solid black line is the best-fit 
re lation obtained using t he weighted x 2 -minimisation procedure 
of Trcmainc et al. I ll2002t) . The dashed line is the best-fit relation 
found using an unweighted least-squares method without taking 
into account error bars. The dotted line is the M — a relation of 
G09. 



avoid systematic errors and for general consistency, we have 
decided to use these results although we are aware that other 
values also exist in the literature. For NGC 1399 and NGC 
5128, G09 give two possible values for the mass. We follow 
their procedure and include both values with a weight of 
1/2 when performing the linear fit, a method also used by 
BT10. For the Milky W ay, we use M BH = 4.3 x 10 6 M Q from 
iGillessen et all (|2009l ) (G09 give M BH = 4.1 x 1O 6 M ). 

For NGC 46 3 6 only two upper limits for A4bh are given 
in iBeifiori et al. I (|20091 ). We do not incorporate these val- 
ues in the determination of the parameters of the relation. 
However, we still include a posteriori NGC 4636 in all of the 
presented figures by taking the mean of the values found in 
Bcifi ori et all |2009) as an upper limit for Mbh. 



3 THE M BH - ct gc RELATION 

In Figure [T] we plot the central black hole mass Mbh ver- 
sus the velocity dispersion ctgc of the GC system with 
the error bars (see Table [1} . To facilitate comparison with 
the previous works of G09, we use the same presenta- 
tion. We assume a relation of the form log (Mbh/Mq) = 
a + j3 log (ctgc/200 km s _1 ) . The parameters of the re- 
lation are calcu l ated by the x 2 - m i n i m isation technique of 
iTremaine et all (l2002h using the IDL MPFITEXY routine 
(Williams, Bur eau, fc Cappellarillioiol i which includes error 
bars in both Mbh and ctgc (weighted fit) and allows the 
determination of the intrinsic scatter eo in Mbh at fixed 
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Table 1. Sample of galaxies with measured line-of-sight velocities of globular clusters and mass of 
the central black hole. 



Galaxy 



Type 



Mbh 

(Mq) 



Ref. 



o"GC 

(km s" 1 ) 



blue 
°GC 

(km s" 1 ) 



rod 
GC 
(km s" 1 ) 



Ref. 



NGC 224 
NGC 524 
NGC 1316 
NGC 1399 



NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 



3031 
3379 
4472 
4486 
4594 
4649 
5128 



NGC 7457 
MW 



M31 



Fornax A 



M81 

M49 
M87 

M104 
M60 

Cen A 



Sb 
SO 
SAB 
El 

Sb 
E0 
El 
El 
Sa 
E2 
SO/E 

SO 
Sbc 



1.5 
8.3 
1.5 
1.3 
5.1 
8.0 
1.2 
1.5 
3.6 
5.7 
2.1 
3.0 
7.0 
4.1 
4.3 



+0.9 
0.3 
+2.7 
1.3 
+0.75 

0.8 
+0.5 

0.66 
+0.7 
0.7 
+2.0 
1.1 
+0.8 

0.58 
+0.2 
0.2 
+ 1.0 
1.0 
+4.4 
4.0 
+0.5 
0.6 
+0.4 
0.2 
+ 1.3 
3.8 
1.2 
1.7 
+0.4 



x 10 8 
x 10 s 
x 10 s 
x 10 9 
x 10 s 
x 10 7 
x 10 8 
x 10 9 
x 10 9 
x 10 s 
x 10 9 
x 10 s 
x 10 7 
x 10 6 
x 10 6 



1 
9 
2 
1 

1 
1 
1 
1 
1 
1 
1 

1 

18 



1M -5 

18613 
20211 
274± 9 



128" 



.+9 
'-9 

312±g 7 
320+i^ 

zi '-16 

150tjj 



hl2 
-12 



129 
197 



+8 
-6 
+39 



29E 



141+15 

342±H 
203 

207±\l 
1491* 



120 



-14 
-14 



121: 

169" 



255 4 



125^ 



265±H 
295lg 
207 

240^0 
;+ 4 



156 



61 



-10 
-10 



10 

11 

12 
13 
14 
16 
17 

3 
15 



REF E RENCES. (1) iGii tekin et al.|l2009l; (2)lNowak et~alll2008l: (3)IChomiuk. Strader. fc Brodie 



Bergond et alj|2006l; (12) ICote et alj|2003t (131 IStrader et aUfipTll; (14). 
Zinnlll996l; (16)lHwang et al. 2008; (17) Woodlcv et al.ll20ld; (18) [Gillcsscn 



2008] ; (4 lShen fc Gebhardl Foi2l ; (5) iBeaslev et alj|2004 (6)lLee et al.ll2008!; (7) TGTudfroor7*et**aT 
200l] ; (8) iRichtler et al.l 12001 (9) jKrainovic et al.l |2009|; 7l0) iNantais fc HuchraT 20ld; (11) 

alj |2003t (13) IStrader et al.ll201ll; (14) iBridges et alj | 2007l; (15) 

et al.ll2009l 



Table 2. Values of the parameter of the relation 
log (Mbh/Mq) = a + /31og (cr/200 km s" 1 ) for the full GC sys- 
tem as well as for the red and blue sub-populations. 







a 









£0 


Full sample 


8.63 


±0.09 


3.76 


± 


0.52 


0.27 


Red GCs 


8.73 


±0.09 


3.84 


± 


0.52 


0.22 


Blue GCs 


8.66 


±0.12 


3.00 


± 


0.72 


0.33 



o-GC ■ The MPFITEXY routine depends on the MPFIT pack- 
age l|Markwardtll2009l ) .We also carry out a standard linear 
least-squares fit without taking into account error bars (un- 
weighted fit) for comparison. We obtain a = 8.63 ±0.09 and 
/3 = 3.76 ± 0.52 for the weighted fit (table 0. For the full 
sample used by G09, these values are a = 8.12 ± 0.08 and 
P = 4.24 ± 0.41. The slope of the relation between Mbh 
and ctgc is consistent with the one obtained by G09. How- 
ever, for the same Mbh the velocity dispersion of the GC 
system is systematically smaller than that obtained for the 
spheroid. The intrinsic scatter we found is £0 = 0.27. As 
a comparison, we have examined the usual Mbh — o- sp h re- 
lation using the galaxies in our sample with values of <j sp h 
taken from G09 and obtained a similar value of eo = 0.27 
with the same fitting procedure. Thus, it seems that, for the 
limited number of galaxies in our sample, Mbh correlates 
equally well with either a sv h or ctgc- 
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Figure 2. Mass A4bh versus velocity dispersion ctqq of the red, 
metal-rich globular cluster system for galaxies in our sample. As in 
Figure [T] circles are data points included in the fitting procedure 
whereas squares are systems which do not contribute to the best- 
fit relation, namely NGC4594 and NGC4636. The lines have the 
same meaning as in Figure [T] 



4 RELATION FOR BLUE AND RED 
GLOBULAR CLUSTER SYSTEMS 

GCs are usually divided into two sub-populations with dif- 
ferent features : the metal-poor GCs which are generally 



older and have a shallower number density profile than the 
metal-rich GCs which consist mostly of younger objects as- 
sociated with the spheroid and lying closer to the centre. 
The two sub-populations are often referred to as the blue 
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Figure 3. Same as Figure [2] but considering the velocity disper- 
sion CTQ^ e of the blue, metal-poor globular cluster system. The 
dash-dotted line is the best-fit obtained by including the Milky 
Way in the fitting procedure. 



and red GCs sub-populations respectively. Here, we examine 
separately the Mbh — cgc relation of these two categories. 

Concerning the red GCs, measurements of ctqq were 
available for 9 systems in our sample (including the red sub- 
population for the Milky Way). The data points as well as 
the best-fit relations are shown in Figure [2] (the same fitting 
techniques detailed in section [3] were used here for the red 
and blue GCs). The relation we find has a = 8.73 ±0.09 and 
ft = 3.84 ±0.52 with an intrinsic scatter £o = 0.22 making it 
a tighter relation than the Mbh — <y sv h for the same sample. 
However, since the Milky Way probes a rather different re- 
gion of parameter space and because its black hole mass is 
known to high accuracy compared to other systems, it could 
potentially introduce a bias in the estimation of the param- 
eters. Nonetheless, we found little difference in the resulting 
parameters even if we removed the Milky Way. We also point 
out that both fitting methods gave almost the same results 
indicating that the parameters are well-constrained. We note 
that our value of ft is close to the value given by G09 for 
early-type galaxies (ft — 3.86) and ellipticals (ft — 3.96). 

For blue GCs (Figure©, we obtain a = 8.66 ±0.12 and 
ft = 3.00 ± 0.72 without including the MW. The intrinsic 
scatter £o = 0.33 is larger than the previous one obtained 
for either the full GC system or for the red GCs and the slope 
is also only marginally consistent with both estimates. Thus, 
it seems that SMBHs and metal-poor GCs are only weakly 
connected. Note that if we include the MW, we obtain a 
higher value for the slope ft = 4.07 ± 0.89 which is more 
consistent with G09 and with the metal-rich GCs but the 
scatter is even larger in this case. 

Finally, we point out that, according to our best-fit re- 
lations for red and blue GCs, there is a preferred value of 
M BH for NGC 1399 and NGC 5128 (Figure [2] and [3} which 



Table 3. Estimated black hole mass using the 
measured velocity dispersion of the red GCs sub- 
population. 



Galaxy 


GC 

(km s" 1 ) 


Ref. 


M BH 

(Mq) 




NGC 1407 


243l?i 


1 


1.12 ±0.42 


x 10 9 


NGC 3923 


200±% 


2 


5.31 ± 2.50 


x 10 8 


NGC 4494 




3 


2.69 ± 2.04 


x 10 7 


M33 


42±s 8 


1 


1.32 ± 1.93 


x 10 6 


LMC 


17 


5 


4.11 ± 5.36 


x 10 4 


REFERENCES. (1) 


iRomanowsky et al 


[200 


(2) iNorris et al.l 12012 


; (3) 


iFoster et al. 


1201 


(4lSchommer et al.lll99ll; (51 lFreemanlll993l 



are Mbh = 3.0 
respectively. 



x 10 8 M and Mbh = 1.3 x 10 9 M, 







5 PREDICTION OF Mbh FOR FIVE GALAXIES 

As we have shown, <Jqq seems to be the best proxy for black 
hole mass for the sample we considered. Thus, we use the 
Mbh — o"gc relation to estimate the black hole mass of galax- 
ies for which crgg is known, namely for NGC 1407, NGC 
3923, NGC 4494, M33 and the LMC (Table [3]). 

• M33 - We obtain Mbh = 1.32 ± 1.93 x lO 6 Af which 
is consiste nt with M33 having no central black hole. We re- 
mark that lGebhardt et alj (|200ll ) have used the Mbh —o- sp h 
relation (a — 8.11 and ft = 3.65) to estim ate the black 
hole m ass in M33. They found 5.6 x 10 4 M g . [Merritt et all 
i|200ll ) find Mbh = 3 x 1O 3 M as an upper limit. They 
conclude that either M33 does not contain a BH or has 
an intermediate mass black hole or that the M — a rela- 
tion cannot be extrapolated to such low masses. Using the 
observed relation bet ween circular vel o city a t large radii 
and black hole mass, Gebhardt et al.l (|200ll ) have calcu- 
lated a mass of Mb h ~ 6 x 10 6 M^. X-ray observations 
l|Foschini et all |2004| ; IZhang et all 120(391 : IWeng et all I2009T ) 
give Mb h ~ lOM^ . Using a lso X-ray data and a different 
method, iDubus et al.l (|2004T ) find a black hole mass in the 
range 10 5 - 10 6 M . Further dat a are needed before a final 
conclusion can be drawn. 

• LMC - The Large Magellanic Cloud, being a satellite 
galaxy, is clearly different from the other systems in the sam- 
ple used to derive the Mbh — o"gc relation. The centre is not 
well defined and its globular clusters are not old. However, 
since the kinematics of the GC system is available (|FreemarJ 
I1993T ). we simply mention that assuming the Mbh — o"gc re- 
lation can be applied it predicts A#bh = 4.1± 5.4 x 10 4 Mq 
which hints towards an intermediate mass black hole but is 
also consistent with the LMC having no central black hole. 



The black hole mass estimates are summarised in Table 



m 



6 CONCLUSION 

In the present paper, we have shown that the velocity disper- 
sion of globular cluster systems projected on the line-of-sight 
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(that is the observed radial velocities) is well-correlated with 
the mass of central black holes, particularly for red (metal 
rich) globular clusters. The slope of the correlation is sim- 
ilar to the one obtained by G09 but the normalisation is 
different, meaning that at fixed black hole mass, the veloc- 
ity dispersion of the globular cluster system is smaller than 
the dispersion of the spheroid. We also show that the rela- 
tion for red clusters has the same slope as that of G09 for 
early-type galaxies. Thus, it seems that the relation with red 
globular clusters is an extension of what happens with o sp h- 
For blue globular clusters, the relation is clearly weaker and 
the slope is also lower than the classical Mbh — a relation. 
Therefore it is not clear whether a direct connection between 
dark matter halos, as probed by blue GCs, and central black 
holes really exists or not. We have also used the tight rela- 
tion with red GCs to estimate the black hole mass in NGC 
1407, NGC 4494, NGC 3923, M33 and the LMC. 
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